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Abstract 

The energy loss effect in nuclear matter , which is another nuclear 
effect apart from the nuclear effect on the parton distribution as in 
deep inelastic scattering process ,can be measured best by the nu- 
clear dependence of the high energy nuclear Drell-Yan process. By 
means of the nuclear parton distribution studied only with lepton 
deep inelastic scattering experimental data, measured Drell-Yan 
production cross sections for 800GeV proton incident on a variety 
of nuclear targets are analyzed within Glauber framework which 
takes into account energy loss of the beam proton. It is shown that 
the theoretical results with considering the energy loss effect are in 
good agreement with the FNAL E866. 
Keywords: Drell-Yan, energy loss,multiple scattering 
PACS:11.80.La,13.85.QK,25.40.-h 

I Introduction 

Propagating of a high energy particle through a nuclear medium is of 
interest in many areas of physics. High energy proton-nucleus scattering has been 
studied for many decades by both the nuclear and particle physics communities W. 
By means of the nuclei, we can study the space-time development of the strong 
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interaction during its early stages, which is inaccessible between individual hadrons. 
The Drell-Yan reaction' 2 ! on nuclear targets provides, in particular, the possibility 
of probing the propagation of projectile through nuclear matter, with the produced 
lepton pair carrying away the desired information on the projectile after it has 
travelled in the nucleus. Only initial-state interactions are important in Drell- 
Yan process since the dimuon in the final state does not interact strongly with 
the partons in the nuclei. This makes Drell-Yan scattering an ideal tool to study 
energy loss' 3 l 

Drell-Yan scattering is closely related to deep-inelastic scattering (DIS) 
of leptons ,but unlike DIS ,it is directly sensitive to antiquark contributions in 
target parton distributions. When DIS on nuclei occurs at x < 0.08, where x is the 
parton momentum fraction, the cross section per nucleon decreases with increasing 
nucleon number A due to shadowing^. Shadowing should also occur in Drell-Yan 
dimuon production at small X2, the momentum fraction of the target parton, and 
theoretical calculations indicate that shadowing in the two reactions has a common 
origin. ^ The energy loss effect is another nuclear effect apart from the nuclear 
effect on the parton distribution as in DIS scattering process. Shadowing and initial 
state energy loss effect are processes that occur in both Drell-Yan reaction and J/ip 
formation. Characterizing the energy loss effect in nuclear matter should further the 
understanding of J/ip production ,which is required if it is to be used as a signal 
for the quark-gluon plasma in relativistic heavy ion collisions. 

In 1999, Fermilab Experiment866(E866) ^ reported the precise measure- 
ment of the ratios of the Drell-Yan cross section per nucleon for an 800GeV proton 
beam incident on Be, Fe and W target at larger values of xi,the momentum frac- 
tion of the beam parton, larger values of xp(~ X\ — rr 2 ),and smaller values of x 2 
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than reached by the previous experiment, Fermilab E772I 7 !. The extended kine- 
matic coverage of E866 significantly increases its sensitivity to energy loss and 
shadowing. 

Recently, M.B.Johnson et al.® and Francois Arleo ^gave theoretical anal- 
ysis of the E866 Drell-Yan experimental data by means of different methods, re- 
spectively. Johnson et al. examined the effect of initial state energy loss on the 
Drell-Yan cross section ratios versus the incident proton's momentum fraction by 
employing a new formulation of the Drell-Yan process in the rest frame of nucleus. 
Francois Arleo carried out a leading-order analysis of E866 Drell-Yan data in nuclei 
according to multiple scattering of a high energy parton traversing a large nucleus 
("cold " QCD matter) studied by Baier et alJ 10 ],in which the multiple soft gluon 
emission from the incoming parton leads to the parton energy loss. 

Since the EMC effect was discovered, various phenomenological mod- 
els have been proposed to investigate the nuclear effect I 11 ' 12 ' 13 - 14 ! .Bickerstaff et 
al.^found that although most of the theoretical models provide good explanations 
for the EMC effect ,they do not give a good description of the nuclear Drell-Yan 
process. Most of them overestimate the Drell-Yan different cross section ratios 
.In previous paper! 16 ' ;We suggested an additional nuclear effect due to the en- 
ergy loss in nuclear Drell-Yan process by means of the Glauber model' 17 ! ,which 
have been extensively employed in nucleus-nucleus reaction with good fit to related 
experiment' 18 ^. It was found that the nuclear Drell-Yan ratio is suppressed signifi- 
cantly as a consequence of continuous energy loss of the projection proton to the 
target nucleon in their successive binary nucleon-nucleon collisions. This suppres- 
sion balances the overestimate of the Drell-Yan ratio only by consideration of the 
nuclear effect from the parton distribution. 
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Recent ly,there were two trials to obtain nuclear parton distributions from 
the existing world experimental data. In 1999,K.J.Eskola et al.(EKS)[ 19 ' suggested 
a set of nuclear parton distributions ,which are studied within a framework of the 
DGLAP evolution .The measurements of F 2 A /F 2 D in deep inelastic I A collisions, 
and Drell-Yan dilepton cross sections measured in pA collisions were used as con- 
straints. The kinematical ranges are 10~ 6 < x < 1 and 2.25GeV 2 < Q 2 < 10 A GeV 2 
for nuclei from deuteron to heavy ones .With the nuclear parton distributions , 
the calculated results agreed very well with the relative EMC and Fermilab E772 
experimental data. In 2001,M.Hirai et al.(HKM)[ 20 l proposed two types of nuclear 
parton distributions which were obtained by quadratic and cubic type analysis. 
,and determined by a x 2 global analysis of existing experimental data on nuclear 
structure functions without including the proton-nucleus Drell-Yan process .The 
kinematical ranges covered 10 -9 < x < 1 and lGeV 2 < Q 2 < !0 5 GeV 2 for nu- 
clei from deuteron to heavy ones . As a result, they obtained reasonable fit to 
the measured experimental data of F 2 . In this report, by means of the EKS and 
HKM nuclear parton distribution functions ,the Drell-Yan production cross section 
rations for 800Gev proton incident on a variety of nuclear targets are analyzed 
by using the Glauber model with taking into account of the energy loss of the 
projective proton through nuclei. 
II Method 

According to Glauber model' 17 ', the projectile proton scattering inelas- 
tically on nucleus (A) makes many collisions with nucleons bound in nuclei. The 

— * 

probability of having n collisions at an impact parameter b can be expressed as 

P(b,n) = n , ( ^ n) , [T(% m ni - Tib)a m ] A -\ (1) 
where o" in (~ 30m6)is non- diffract ive cross section for inelastic nucleon-nucleon col- 
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lision,and T(b) is the thickness function of impact parameter b. For collisions of 
nucleons which are not polarized, the collisions does not depend on the orientation 
of b , and T(b) depends only on the magnitude |6| = b. So, we could consider only 

— * 

this case of T(b) = T{b). In an nucleon-nucleus collision without impact param- 
eter selection, the number of nucleon-nucleon collisions n (for n — 1 to A) has a 
probability distribution 

P(„) = j SP ^ . (2) 
E fdbP{n,b) 

n=l 

In following calculation , the thickness function can be conveniently written as' 18 ' 

- f ^rexp(-b 2 /2p 2 A ), A<32, 
T{b) = , ^a^_ HA! ^ (3) 
\ ^rJR 2 A-V6(R A -\b\), A>32. 

Here Ra = tqA 1 ^ 3 is the radius of a colliding nucleus with r = 1.2 fm, 9 is a step 
function , and f3 A = 0.606A 1 / 3 . 

In multiple-collision Glauber model, the basic process is nucleon-nucleon 
collision for proton-nucleus Drell-Yan process. The leading-order contribution to the 
Drell-Yan process is quark-antiquark annihilation into a lepton pair. The annihila- 
tion cross section can be obtained from the e + e~ — > fi + fi~ cross section by including 
the color factor | with the charge e 2 for the quark of flavor /. 

= e 2 f 5(s-M 2 ), (4) 

dM 9M f y h y ' 

where s = xix 2 s,is the center of mass system (c.m.system) energy of qq collision 
,a;i(resp.x 2 )is the momentum fraction carried by the projectile (resp.target)parton, 
y/s is the center of mass energy of the hadronic collision, and M is the invariant 
mass of the produced dimuon. The hadronic Drell-Yan differential cross section is 
then obtained from the convolution of the above partonic cross section with the 



quark distributions in the beam and in the target hardron: 

' /<T =^^^E^[^i)^)+^(xi)^(x 2 )], (5) 



dxidM ' 9M xis , 

where K is the high-order QCD correction, a is the fine-structure constant, and 
(x) and (f^ (x) are the quark and anti-quark distributions in the proton (nu- 
cleon in the nucleus A) .In addition ,one has the kinematic relations , 

M 2 

XiX 2 = ,X F =X!-X2, (6) 

s 

with the Feynman scaling variable 

x F = —=, (7) 
V s 

where pi is the longitudinal momentum of the virtual photon , 

Now let us take into account of the energy loss of the projectile proton 
moving through nuclei. In proton-nucleus Drell-Yan collision, the incident proton 
passes through the nucleus befor producing the high Q 2 dimuon pair. On the one 
hand ,the projectile proton may interact spectator nucleon bound in nuclei ,in which 
soft (nonperturbative) minimum bias collisions may occur before making the final 
lepton pair. As a result, the projectile proton imparts energy to the struck nucleon 
and therefore must loose energy. On the other hand, the projectile proton, which 
travels through nucleus with existence of multiple-collisions , may emit soft gluon 
between two bias collisions in the nucleus ,and must consequently experience an 
energy loss. Thus the energy loss from two above aspects in multiple collisions can 
induce the decrease of cm. system energy of the nucleon-nucleon collision produc- 
ing dimuon ,and affect the measured Drell-Yan cross section. After the projectile 
proton has n collisions with nucleons in nuclei, suppose for conveniently calculation 
that the c.m.system energy of the nucleon-nucleon collision producing dimuon can 



be expressed as 

>/s 7 =< s /i-(ra-l)A<v/i, (8) 

where A^/s is the c.m.system energy loss per collision in the initial state . There- 
fore, the cross section for Drell-Yan process can be rewritten as 

Here the rescaled quantities are defined as 

x' 12 = r s x 1>2 , (10) 

with the c.m.system energy ratio: 
because of 

/ 2 Pi / 10 N 

x F = -= = r s x F . (12) 



It is obvious that r s is always greater than one. Combining above ingredients, The 
measured Drell-Yan cross section in proton-nucleus collision experiments can be 
expressed as 

III Results and Discussion 

In order to compare with the experimental data from E866 collaboration' 6 l,we 
introduce the nuclear Drell-Yan ratios as: 

Ra " aAxi) = i^m-y (14) 

The integral range on M is determined according to the E866 experimental kine- 
matic region. In our theoretical analysis ,x 2 is calculated with the Drell-Yan differ- 



ential cross section rations Ra 1 /a 2 as 

( TDdata __ ptheo \2 
v 2 _ \ - v A1/A2J n A 1 /A 2 ,j) n ,x 
A ^ /perr \2 ' V ' 

where the experimental error is given by systematic errors as R e X[/A 2 j> anc ^ ^aI^a^ j( 
R t M°A 2 j) indicates the experimental data (theoretical values ) for the ratioi?^/^- 

If the EKS^ 19 ! nuclear parton distribution functions are used together with 
CTEQ (the coordinated Theoretical Experimental Project on QCD )' 21 ^ parton 
distribution functions in proton, obtained x 2 value is x 2 — 51.4 without energy 
loss effects for the 56 total data points. The x 2 P er degrees of freedom is given 
by X 2 '/d.o.f. = 0.918. It is shown that theoretical results agree very well with the 
E866 experimental data ,which results from EKS parametrization of nuclear parton 
distributions studied with including the Drell-Yan process. 

In addition, we consider also using HKM^ 20 ! nuclear parton distribution 
functions together with MRST t 22 ' parton distribution functions in proton. The cal- 
culated results with HKM cubic type of nuclear parton distribution are shown 
in Fig.l and Fig.2. which is the Drell-Yan cross section ratios for Fe to Be 
and W to Be as functions of x\ for various interval of M, respectively. The solid 
curves are the ratios with only the nuclear effect on the parton distribution as 
in DIS scattering process, and the dotted curves correspond to an energy loss ef- 
fect: Ay/s = 0.18GeV with nuclear effect on structure function . Obtained x 2 
values are x 2 — 56.39 with energy loss effect (i),x 2 = 143.74(m) without energy 
loss effect, respectlyfor the 56 total data points. The x 2 P er degrees of freedom 
are x 2 /d-o.f. = 1.007(i), 2. 567(H), respectively. Employing the HKM quadratic 
type of nuclear parton distribution ,x 2 values are x 2 — 56.81 with energy loss 
effect, x 2 = 143.88 without energy loss effect, respectly.The x 2 P er degrees of free- 
dom are x 2 /d-o.f. = 1.015,2.569 , respectively. This implies that the observed 
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energy loss of the incident proton is At/s = 0.18GeV at 1<t, which is almost equal 
to our previous resultsAy's = 0.2GeVt 16 l From comparison with the experimental 
data , it is found that our theoretical results with energy loss effect are in good 
agreement with the Fermilab E866. 

In summary, we have made a leading-order analysis of E866 data in nuclei 
within the framework of the Glauber model by taking into account of the energy 
loss effect of the beam proton. In continuous multiple collisions ,the energy loss ef- 
fect in initial state causes the suppression of the proton-induced nuclear Drell-Yan 
cross section . It is found that the theoretical results with energy loss are in good 
agreement with the Fermilab E866 experiment by means of the parametrization 
of nuclear parton distributions studied without nuclear Drell-Yan process. This is 
another nuclear effect apart from the nuclear effect on the parton distribution as in 
DIS scattering process. The nuclear effect on structure functions and initial state 
energy loss can occur in both Drell-Yan production and J/tp formation. Hence, 
these researches should also further the understanding of J/ip production. Con- 
sidering the existence of energy loss effect in Drell-Yan lepton pairs production, we 
think that the determination of nuclear parton distribution functions should not 
include Drell-Yan experimental data.Although there are abundant data on electron 
and moun deep inelastic scattering, valence quark distributions in the small x re- 
gion and the anti-quark distributions are difficultly determined .At this stage ,only 
valence quark distributions in medium x region can be relatively well determined. 
It is well considered that the precise nuclear parton distributions must be known 
in order to calculate cross sections of high energy nuclear reactions accurately and 
find a signature of quark-gluon plasma in high energy heavy-ion reactions. We sug- 
gest precise neutrino scattering experiments ,which can provide a good method 
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for measuring the F 2 (x, Q 2 ) and xF 3 (x,Q 2 ) structure functions. Using the average 
of xF£ A (x, Q 2 )and xF£ A (x, Q 2 ) , the valence quark distribution functions can be 
well determined. Combining the lepton inelastic scattering data with the neutrino 
scattering experiments, valence quark and anti-quark distribution functions will be 
obtained in the future, which makes us good understanding the energy loss effect 
in high energy nuclear Drell-Yan collisions. 

Acknowledgement: This work is partially supported by Major State Basic 
Research Development program (under Contract No. G20000774), by the Na- 
tional Natural Science Foundation of China(No.l9875024,19835010,10175074), and 
by Natural Science Foundation of Hebei Province(No. 103143) 

References 

[1] W.Busza R Ledoux, Ann.Rev.Nucl.Part.Sci.,38(1989)119. 

[2] S.Drell and T.M.Yan, Phys.Rev,Lett. ,25(1970)316. 

[3] G.T.Garvey and J.C.Peng ,Phys.Rev.Lett., 90(2003)092302. 

[4] M.Arneodo.et.al.(EMC), Nucl.Phys,B441 (1995)3. 

[5] S. J. Brodsky,A.Hebecker,E.Quark.Phys.Rev.D. 55. (1997)2584. 

[6] M.A.Vasiliev,et.al.(E866),Phys.Rev.Lett. 83(1999)2304. 

[7] D.M.Adle et al.(E772),Phys.Rev.Lett., 64(1990)2479 

[8] M.B. Johnson, B.Z.Kopeliovich and 
I.K.Potashnikova,Phys.Rev.Lett. ,86(2001)4483 

[9] F.Arleo,Phys.Lett.,B532(2002)231 

10 



[10] R.Baier,Yu.L.Dokshitzer,A.H. Mueller, S.Peigne, and 
D.Schiff,Nucl.Phys.,B484(1997)265 

[11] C.A.Garacia C anal, E.M.S ant angle, and 

H.Vucetich,Phys.Rev.Lett., 53(1984)1430 

[12] F.E.Close,R.L.Jaffe,R.G.Roberts and G.G.Ross,Phys.Rev.,D31(1985)1004 

[13] G.L.Li,K.F.Liu,G.E.Brown,Phys.Lett.,B213(1998)531 

[14] Zhenmin He,Xiaoxia Yao,Chungui Duan,Guanglie Li,Eur.Phys. J.C4(1998)301 

[15] R.P.Bickerstaff,M.C.Birse,G.A.Miller,Phys.Rev.,D33(1986)322 

[16] Jianjun Yang and Guanglie Li,Eur.Phys.J.C5(1998)719 

[17] R.J. Glauber, Lecture in Theoretical Physics, edited by W.E.Brittin 
and L.G.Dunham(New.York,1959),vol.l,p.315. 

[18] C.Y.Wong,Phys.Rev.D30(1984)961;C.Y.Wong,Introduction to 
High-Energy Heavy-Ion Collisions, 

(World Scientific, Pulishing, Coe Singa Pte.Ltdl.,1994),p.249. 

[19] K.J.Eskola,V.J.Kolinen and C.A.Salgado(EKS),Eur.Phys.J.C9(1999)61. 

[20] M.Hirai,S.Kumano,M.Miyama(HKM) ,Phys.Rev.D64(2001)034003. 

[21] H.L.Lai, et al.(CTEQ),Eur.Phys.J.C5(1998)461. 

[22] A.D.Martin, R.G.Roberts, W.J.Stirling ,and R.S.Thorne(MRST), 
Eur.Phys.J.C4(1998)463. 

Figure caption 
11 



Fig.l The nuclear Drell-Yan cross section ratios Ra 1 /a 2 ( x i) 011 F e to Be for various 
intervals M. Solid curves correspond to nuclear effect on structure function . Dotted 
curves show the combination of shadowing and energy loss effect ( Ay/s = 0.18GeV) 
with HKM cubic type of nuclear parton distributions .The experimental data are 
taken from the E866[6]. 

Fig. 2 The nuclear Drell-Yan cross section ratios Ra 1 /a 2 ( x i) 011 W to Be for 
various intervals M. The comments are the same as Fig.l 
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